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We exhibit a resonance mechanism of amplification of density perturbations in inflationary mo- 
dels, using a minimal set of ingredients (an effective cosmological constant, a scalar field minimally 
coupled to the gravitational field and matter), common to most models in the literature of inflation. 
This mechanism is based on the structure of homoclinic cylinders, emanating from an unstable 
periodic orbit in the neighborhood of a saddle-center critical point, present in the phase space 
of the model. The cylindrical structure induces oscillatory motions of the scales of the universe 
whenever the orbit visits the neighborhood of the saddle-center, before the universe enters a period 
of exponential expansion. The oscillations of the scale functions produce, by a resonance mechanism, 
the amplification of a selected wave number spectrum of density perturbations, and can explain the 
hierarchy of scales observed in the actual universe. The transversal crossings of the homoclinic 
cylinders induce chaos in the dynamics of the model, a fact intimately connected to the resonance 
mechanism occuring immediately before the exit to inflation. 
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The existence of an inflationary phase in the early stages of our Universe has become one of the paradigms of 
modern cosmology, and is now being subject to experimental verification through the crucial measurements of small 
scale anisotropies in the cosmic background radiation. The basic physical ingredient for this inflationary phase is the 
existence of a scalar field - the inflaton field, the vacuum energy of which plays the role of a cosmological constant, 
engendering via the gravitational dynamics an exponential expansion in the comoving scales of the universe. This 
phase is considered to have evolved from a pre-inflationary phase just exiting the Planck era, the matter content of 
which may be thought to be modeled by radiation plus a scalar field minimally coupled to the gravitational field. 
To be considered successful, the inflationary paradigm must accomodate models in which initial homogeneities are 
generated, leading to the formation of the structure present in our actual universe. In the literature of inflation, it 
has been claimed that small fluctuations in the scalar field may give rise to an almost scale-free spectrum of adiabatic 
density perturbations (the Harrison-Zeldovich spectrum), which is the standard spectrum used in explaining galaxy 
formation. However the number of distinct large scale objects and structures in our universe, as galaxies, cluster of 
galaxies, quasars, filaments and walls, is very large and the sizes of these objects form a hierarchy of scales that is not 
described by the flat spectrum of the initial perturbations Q. In other words, the scale- free adiabatic perturbations 
appears not to be sufficient to produce the large-scale structure of the observed universe. We are therefore led to 
pose the following basic question: would it be possible, in the realm of the inflationary paradigm, to obtain adiabatic 
perturbations with a spectrum having several maxima, which would help to explain the origin of the hierarchy of 
scales? 

Our purpose in this essay is to give a positive answer to the above question, and show that the ingredients present 
in most inflationary models may engender a resonance mechanism, in the beginning of, or immediately before the 
inflationary period, that amplifies scalar field perturbations with a selected wave number spectrum. We emphasize the 
simplicity of this mechanism, which deals with a minimum set of ingredients only and does not demand the introduction 
of further fields, fine tuning of parameters and special initial conditions. Even if the universe inflates afterwards, the 
relative rate of amplitudes produced after the resonance amplification would be maintained as an imprint in the initial 
spectrum of density fluctuations. We will also show that this resonance amplification mechanism has a fundamental 
connection to the chaotic behaviour of the model, and its chaotic exit to inflation 0] . Similar resonance mechanisms 
have been previously applied for particle production during reheating in the realm of inflationary models J3|. 

Let us consider the following configuration, that may constitute a simple description of a preinflationary phase 
of the universe in most of the inflationary models discussed in the literature. We start from a Bianchi IX model 
characterized by two scale factors A(t) and B(t), with a perfect fluid and a cosmological constant A playing the role 
of the vacuum energy of the scalar field. For simplicity, we consider the case of dust, but the basic features are not 
altered for other types of perfect fluid like radiation, for instance. According with Ref. 0, the full dynamics of the 
models is governed by the Hamiltonian 

H(A,B,P A ,P B ) = ^^-^ + 2A-^-2AAB 2 -Eo = 0, (1) 

where Pa and Pb are the momenta canonically conjugated to A and B, respectively, and Eq is a constant proportional 
to the total matter energy of the models. The associated dynamical system has one critical point E of saddle-center 
type in the finite region of the phase space. The phase space has two critical points at infinity, corresponding to 
the stable and unstable de Sitter solutions, and the scale factors A(t) and B{t) approach the de Sitter attractor 

exponentially, A = B = eV A / 3t . The dynamical system generated by (1) admits the invariant manifold M. (A = B, 
Pa = ^r) in which all orbits are Friedmann models with cosmological constant and dust. The critical point E is 
contained in VW, and that the separatrices S constitute boundaries between isotropic models that collapse or escape 
to the de Sitter configuration. The models to be considered in this essay are the ones restricted to a neighborhood 
of the invariant manifold M and can be understood as describing the early stages of inflation in which anisotropy is 
present, in the form of small fluctuations. 

The basic point of the amplification mechanism to be discussed here is the structure of homoclinc cylinders present 
in the neighborhood of the saddle-center. To see this we make use of a theorem by Moser Q] which states that it 
is always possible to find a set of canonical variables, say (<?]., <?2 1 Pi jP2)j such that, in a small neighborhood of a 
saddle-center at the origin, the Hamiltonian is separable as 

VA, a 



H(q u q 2 ,p 1 ,p 2 ) = —(pl- eg) - ^2K{p{ + qf) - E + E cr a 0, (2) 

where E cr is the energy associated to the saddle-center, with the two approximatte constants of motion given by the 
partial energies 
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E hyp = ^{p\- q 2 2 ), E rot = V2 A (p( + q{). (3) 
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The energies Eh yp and E rot will be referred to as the hyperbolic motion energy and rotational motion energy of the 
system about E, respectively. The rotational motion corresponds to periodic orbits in a linear neighborhood of E. For 
zero energy of the hyperbolic motion, we have the linear stable V s and unstable V u one-dimensional manifolds. The 
separatrices S in the invariant manifol M. are actually the non-linear extension of V s and V u . The direct product of a 
periodic orbit with V s and V u generates, in the linear neighborhood of E, the structure of stable and unstable cylinders, 
which coalesce into the periodic orbit for large positive and negative times, respectively ||. A general orbit which 
visits the neighborhood of E has an oscillatory approach to the cylinders, the closer as the energy of the hyperbolic 
motion goes to zero. Away from the linear neighborhood of E, the linear approximation is no longer valid. Higher 
order terms become important for the dynamics, and the non-integrability of the system results in the distortion 
and twisting of the cylinders. The stable cylinder and the unstable one will cross each other transversally, producing 
chaotic sets || in the neighborhood of the invariant manifold M. ; orbits with initial conditions taken on these chaotic 
sets will be highly sensitive to small perturbations of the initial conditions. The chaotic character manifests itself 
physically such that a small perturbation of initial conditions taken on these sets would change an orbit from collapse 
to escape into the de Sitter phase, and vice-versa. Also, it is always possible to select initial conditions in these sets 
such that the scale factors oscillate an arbitrary fixed time before escaping to a de Sitter phase. 

The physical effects related to the growing of inhomogeneous perturbations occur during the oscillatory phase when 
a given orbit approach E before to collapse or to escape. We consider a given orbit initially in a very small sphere of 
initial conditions about the invariant manifold M.. In this regime the scale factors oscillate in a small region about 
the critical point E. To avoid technical dificulties concerning the perturbative formalism in anisotropic backgrounds, 
and motivated by the fact that in the overall dynamics of perturbations these small anisotropy fluctuations can be 
averaged, we introduce the average scale factor l(t) of the universe defined by 

1 = {AB 2 ) 1/3 - (4) 

This approximately corresponds to a closed Friedamnn-Robertson- Walker backgound with l(t) as the scale factor, 
either in the oscillatory regime and in the inflationary phase. In a very small neighborhood of E the approximate 
analytical solution for the average scale factor is 

l(t) sa l + Ci e t/2l ° + C 2 e~ t/2l ° +ecos (2wt) (5) 

where Iq — l/v4A, u> = \/2 A, and e, C\ and Ci are very small arbitrary constants of the same order of R. The 
exponential terms are associated to the hyperbolic energy mode, while the last term is associated to the rotational 
energy mode. The pure oscillatory phase is characterized by G\ = C2 = (on the cylinder), which is a very special 
class of solution, hardly realized by physically reasonable universes. The escape to inflation and the collapse depend 
on the relative signs of C\ and C2, but the smaller are these constants the larger is the time that the orbit spends 
in the oscillatory phase. Nevertheless it must be stressed that in this oscillating phase the universe has periods 
of accelerated expansion followed by accelerated contraction, and the escape to inflation actually corresponds to a 
progressive increase of the inflationary phase. 

The next step is to study the behavior of inhomogeneous perturbations of matter d p and scalar field 5 ip in this 
background. Before presenting the evolution equations for the perturbed quantities, it is important to stress that 
Sip will give rise to second order perturbations of the energy of the scalar field, 5 p v , since ip^ » V'(ip^) ~ 0, 
with ip 1 -^ m constant being the scalar field of the background corresponding to an extremum of the potential. As a 
consequence, there will be no coupling between S p and 5 p Vl once the former represents first order deviations. 

Introducing the contrast density X = ^f~i one can show that its evolution equation is given by 

i 1 

Xk + 2jx k - -pxk = 0- (6) 

Here Xk is the mode associated to the scalar harmonic basis Q(fe)(x) which satisfies h^ u Q(fc)(x) iAi; „ = ^ Q( /c )(x), h^ u 
being the projector into the spatial surface, and k 2 — n 2 — 1 with n > 1 assuming integers values, p = E /l(t) 3 is the 
energy density of the background. In the oscillatory phase about E (Eq. (||) with C\ = C2 = 0), the above equation, 
after a proper rescaling of the variable Xk, assumes the form 

Xk + (-q + a cos (2wt)) Xk = 0, (7) 

where q is a positive constant, and \a\ << q. This is the Mathieu equation inside the parametric region of instability, 
in which all modes k are equally amplified. If the universe, after this oscillating phase, undergoes to a successful 
inflation, I ps exp y/A/3 1, and the equation for all modes Xk reduces to 
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Xk 



0. 



(8) 



As a consequence all modes Xk will decrease asymptotically, or in other words, all perturbations of matter produced 
during the oscillatory phase will be washed out in the inflationary regime, and cannot therefore be considered as seeds 
of structure formation. 

Now we consider the evolution of perturbations of the scalar field <p(x,t) — <p^(t) + Qr k \(x) 5 <^fe(t), where the 
background ip^(t) m constant. The equation of motion for the perturbed modes of the scalar field is 



(8<p k ) +3 l - {8 Vk ) + (V'V 0) ) + S<p k = 0. 



(9) 



Introducing the variable 8 k 
regime, it follows 



ft/ 2 5 ip k into the above equation as well as the approximate solution in the oscillatory 
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~fT + ~1T- cos(2wt 



5 k = 0, 



(10) 



that is a typical Mathieu equation. As an important property of solutions of this equation is the existence of an 
instability zone in the parameter space in which the solutions grow exponentially as 8kexp{vkt). This instability is 
known by parametric resonance. For small e, the condition for parametric resonance is given by 
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f + A(e) 



(11) 



where j is an arbitrary integer and A(e) is a small adimensional number dependent of e. This means that only the 
modes with k satisfying the above relation grow exponentially. These modes are determined by the ratio between 

the mass of the scalar field V"((p^), and the vacuum energy of the scalar field. The ratio — ^ — - is a fundamental 
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parameter in inflationary models; since the resonance mechanism is believed to be general, the order of this parameter 
could be evaluated through Eq. ( |ll| ) once the spectrum is experimentally determined. It will be quite ilustrative to 
choose standard values used for the afore mentioned constants as V"((f( ') ~ (10~ 12 — 10~ 14 ) rr? { and A = V{(p^) ~ 
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Pi 



where m p i is the Planck mass. We obtain a very small number of order of 10 
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l ) rn 2 ; and A 
10~ 14 that can be negleted. 



Thereupon, during the oscillatory phase a selected spectrum of modes, centered about the integers j = 2 k , grow 
rapidly, whereas the remaining modes are maintained stable. Suppose that, after this regime, the universe escapes to 
inflation. The equation of motion (^J) assumes its usual form as discussed in several works. As we know the perturbed 
modes that cross the horizon are not affected by microphysical processes, being in this way frozen until they reenter 
the horizon after inflation. However, due to the mechanism of parametric resonance, the modes which cross the 
horizon for the first time have not the same amplitude, and are maintained during the kinematical evolution until 
they reenter the horizon. At the end, the resulting spectrum of adiabatic perturbations is not a flat one, in the sense 
that there are maxima dictated by Eq. (11). The non-flat spectrum of inhomogencous fluctuations produced by the 
above mechanism can explain the hierarchy of scales (galaxies, clusters, etc) as observed in the actual universe. 
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